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Abstract In the "braneworld scenario" ordinary standard model matter 
and non-gravitational fields are confined by some trapping mechanism to the 
4-dimensional universe constituting the D3-branes which are embedded in a 
(4 -|- n)-dimensional manifold referred to as the 'bulk' (n being the number 
of extra dimensions). The notion of particle confinement is necessary for 
theories with non-compact extra dimensions, otherwise, the particles would 
escape from our 4-dimensional world along unseen directions. In this paper, 
we have considered a five-dimensional warped product space-time having 
an exponential warping function which depends both on time as well as on 
the extra coordinates and a non-compact fifth dimension. Assuming that 
the lapse function may either be a constant or a function of both time and 
of the extra coordinates, we have studied the nature of the geodesies of test 
particles and photons and have analyzed the conditions of stability in this 
geometrical framework. We have also discussed the possible cosmology of 
the corresponding (3 -I- l)-dimensional hypersurfaces. 

Key words Braneworlds • Geodesic motion • Cosmology. 



1 Introduction 



According to the string theory postulate, standard model matter and non- 
gravitational fields are confined by some trapping mechanism to the 4- 
dimensional universe constituting the D3-branes (4-dimensional timelike hy- 
persurfaces) that are embedded in a (44-n)-dimensional manifold referred 
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to as the 'bulk' (n being the number of extra dimensions). This has led 
to a renewed interest in extra-dimensional theories of gravity. The success 
of the Kaluza-Klein [1] theory in particle physics led several workers to 
construct numerous models of non-compact higher-dimensional theories of 
physics [2] -[6], and eventually to the so-called "braneworld scenario" [7]- 
[9] . In the braneworld models with non-compact extra dimension, particles 
and non-gravitational fields are assumed to be confined to the branes. At 
low energies, gravity is localized at the brane along with the particles but 
at high energies gravity "leaks" into the higher-dimensional bulk and is 
propagated therein. These "braneworld" models have been used to address 
several issues, one of which is to explain why the observable universe is 
found to be 4-dimensional [10], [11]. Consequently, it is necessary to deter- 
mine how closely the corresponding apparent 4-dimensional world resemble 
the observed world [12], and thus we need to verify the extent to which 
these theories satisfy the geodesic postulates concerning the motion of test 
particles and light rays. 

The notion of particle confinement is necessary for theories with non- 
compact extra dimensions, as otherwise, the particles would escape from 
our 4-dimensional world along unseen directions. In the classical context, 
confinement of a test particle to the brane eliminates the effects of extra di- 
mensions, thereby rendering them undetectable. In the braneworld scenario, 
the stability of the confinement of matter fields at the quantum level is made 
possible by assuming an interaction of matter with a scalar field. It appears 
that non-gravitational forces acting in the bulk and orthogonal to the space- 
time are necessary to keep the test particles moving on the brane, the source 
of these confining forces being interpreted in different ways [13]-[15]. The 
alternative explanation in case of higher-dimensional theories involve the 
invocation of geometrical mechanisms to support lower-dimensional con- 
finement. In this case, confinement is due purely to the classical gravi- 
tational effects, without requiring the presence of brane-type confinement 
mechanisms. This has led to the investigation of warped product spaces and 
their geometrical properties [16]. Randall and Sundrum [10], [11], achieved 
it through the construction of five-dimensional warped product spaces using 
an exponential warp factor in a non-factorizable metric furnished with mir- 
ror symmetry, even when the fifth dimension was infinite. In case of warped 
product spaces, it has been observed that a general qualitative analysis of 
the behavior of massive particles and photons in the fifth dimension can 
be made from the knowledge of the warping function. A general pictiire of 
these geodesic motions has been obtained by using the natural decoupling 
occurring between motions in the brane and motion in the fifth dimension 
in case of such spaces [17], [18]. For a thin asymmetric braneworld, bulk 
and brane geodesies does not coincide in general. However, if appropriate 
energy conditions are satisfied by the matter confined to the brane, then 
in presence of mirror-symmetry, test particles can be confined by gravity 
to a small region about the brane [19], [20]. The Z2-symmetric braneworlds 
have their apparent and bulk geodesies coinciding on the brane. Although 
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null bulk geodesic motion in the RS2-type brancworlds, as well as in the 
static universe in the bulk of a charged topological AdS black hole have 
been studied [21], the same is not found for timelike geodesies. 

Another important program in the study of higher dimensional models 
is the cosmological interpretation of the corresponding 4-dimensional ge- 
ometry. By use of the Einstein equations, it has been shown that [22] the 
embedding of a surface in a flat space of co-dimension one imposes the re- 
striction that the surface has a constant curvature, if its dimension is n > 2. 
The effective equations for gravity in four dimensions were obtained by Shi- 
romizu et al [23] and subsequently cosmological solutions have been studied 
by some authors [24], [25]. 

In this paper, we have considered RS-type braneworlds with the bulk 
in the form of a five- dimensional warped product space-time, having an ex- 
ponential warping function which depends both on time as well as on the 
extra coordinates and a non-compact fifth dimension. We know that the 
exponential warp factor reflects the confining role of the bulk cosmological 
constant [26] to localize gravity at the brane through the curvature of the 
bulk. It is, therefore, possible that such localization may include some time- 
dependence and hence the choice. Employing the technique used by Dahia 
[17], we have obtained a mechanism for the confinement of geodesies on 
co-dimension one hypersurfaces, considering that the confinement is purely 
due to the classical gravitational effects. Wo have boon able to obtain the 
description of the geodesic motions by using the natural splitting occurring 
between the motion in the extra dimension and the motion in the four- 
dimensional hypersurfaces. Such splitting helps us to use the phase space 
analysis to determine the nature of the geodesic motions in the neighbour- 
hood of the hypersurfaces and analyze the conditions of stability. Further, 
we have assumed the lapse function to be either a constant or a function 
of both time and of the extra coordinates. This is possible since, in our 
case, the square of the lapse function is identical to the metric coefficient 
for the fifth dimension, which may therefore be of such a type. Finally, 
the cosmological interpretations of the corresponding (3 + l)-dimensional 
hypersurfaces have also been discussed on the basis of the field equations. 

2 Geometric Construction of Five-Dimensional Warped Product 
Spaces 

A warped product space [27], [28], is constructed as follows: Let us consider 
two manifolds (Riemannian or semi-Riemannian) (M™,ft,) and {M",h) of 
dimensions m and n, with metrics h and h respectively. Given a smooth 
function / : M" — > 5R (henceforth called the warping function), we can 
build a new Riemannian (or semi-Riemannian) manifold (M, g) by setting 
M = X M", which is defined by the metric g = e^^h^h. Here (M, 
g) is called a warped product manifold. The case dim M = 4 corresponds 
to (M, g) being a spacetime, and is called a warped product spacetime (or 
simply warped spacetime). 
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In local coordinates z^, the line element corresponding to the bulk metric 
is denoted by 

dS^ = QABdz'^dz^. (1) 

The class of warped geometries considered by us is represented in general 
by the line element of the form 

dS'^ = e^fK0dx°'dx^ + Kbdy^dy^, (2) 

where the extra coordinates are represented by y", the coordinates on the 
m-dimensional submanifold by x'', / is a scalar function, = ha0{x) is 
the warp metric on the submanifold of dimension m and hab is the metric 
representing the extra dimensional part. 

In this paper, we consider m = 4 and n = 1, so that M = x M^, 

where M'' is a Lorentz manifold with signature (H ). The canonical 

metric due to Mashhoon et. al. [29], [30] is more general than the Randall- 
Sundrum metric and encompasses all the metrics usually considered in the 
braneworld and Induced Matter Theory (IMT) approaches [6],[31]-[34]. For 
this, the physical metric in 4-dimension is assumed to be conformally related 
to the induced metric [35], i.e. 

dS'^ = Qds^ + e^'^dy'^, (3) 

where is the "warp" factor that satisfies the condition 12 > and the 
scalar <P which normalizes the vector normal to the hypersurfaces y = 
constant is known as the lapse function. We assume that the branc is de- 
fined by the y = constant hypersurface, where y is a Gaussian normal 
coordinate orthogonal to the brane, representing the fifth dimension, which 
is non-compact and curved (warped) [26], [11]. 

In general, the warp factor is assumed to be a function of extra coordi- 
nates only, but here we assume that it is a function of both time, as well as 
of the extra coordinates. As in the RS models, gravity is localized on the 
brane through the curvature of the bulk. We know tliat the bulk cosmologi- 
cal constant acts to squeeze the gravitational field closer to the brane, with 
the exponential warp factor reflecting the confining role of the bulk cos- 
mological constant [26] . Mathematically, the time dependence of the warp 
factor docs not affect its smooth nature. Physically, it takes into account 
the possibility that the confining role of the bulk cosmological constant and 
the curvature of the bulk may have some dependence on time, in addition 
to their dependence on the extra dimensional coordinate. When the extra 
dimension is spacelike, we have e = — 1. In our case, = p = Qyy (which 
means that the square of the lapse function is identical to the metric coef- 
ficient for the fifth dimension) and 1? = e^^^*''^^ . The expression for the line 
element therefore takes the form 



^52 ^ e'^f^^'^^Kpdx'^dx^ - pdy'^. (4) 
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Here, p may either be a constant or a function of coordinates. We shall 
consider two cases, viz. p = 1 and p = p{t, y) (for example, the second type 
can be found in [36]). To take into account the time-dependence of the warp 
factor, we assume the scalar function /(t, y) to be of the form 

f{t,y) = at + l\n{cosh{cy)) (5) 

where a, I and c are constants. We will not make any assumption on the 
topology of the extra dimension, so that, in principle — oo < y < +oc. 



3 Geodesic Motion 

Dahia [17] demonstrated oscillatory confinement of massive particles and 
light rays for branes of finite thickness in the context of 4 + 1-dimensional 
warped product spaces, whereby, the particles in the four-dimensional hy- 
persurface can oscillate about the hypersurface, while remaining close to it. 
They have shown that such behaviour can occur for large classes of bulks 
which possess warped product geometries. Writing the geodesic equations 
for warped product spaces, they showed that the equation that describes the 
motion in 5D decouples from the rest. They made a qualitative analysis of 
the motions by rewriting the geodesic equation in the fifth dimension as an 
autonomous planar dynamical system and then employed the phase plane 
analysis to study the motion of particles with nonzero rest mass and pho- 
tons respectively. Consequently they could draw very general conclusions 
about the possible existence of confined motions and their stability in the 
neighbourhood of the hypersurfaces. In our analysis, which now follows, we 
have used this method to arrive at the results. 

The equations of geodesies in the 5-dimensional space M is given by 

where A is an afHne parameter and ^^^F^q are the 5-dimensional Christoffel 
symbols of the second kind defined by ^^^T^q = \g^^{gDB,c + 9dc,b — 
9bc,d)- If we denote the fifth coordinate by y and the remaining "space- 
time" coordinates by x^^ i.e. z^ = {x^,y), it can be shown that the 
4-dimensional part of the geodesic equations (6) can be rewritten in the 
form [17] 

+ -'a/3 j\ j\ - ? ' \') 



rfA2 dX dX 



where 



(8) 

with (^^r^s = y^''{gua,0 + 9wfS,a - 9a0,v)- 
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We consider that the bulk spacetime is a 5-dimensional warped product 
spacetime in the braneworld scenario, represented by the line element (4) 
with the warping function given by (5). Let us now assume that this bulk 
space-time is foliated by a family of hypersurfaces defined by the equation 
y = constant. The geometry of each such leaves of foliation denoted by, say 
y = yo will be determined by the induced metric 



ds'^ = ga0{x,yo)dx"dx^ = e'^^^*'y°^K0{x)dx'^dx^ = r)c,pdx°'dx^ . (9) 

The extrinsic curvature of such a hypersurface is given by 2Kap = 
— (^g^). In that case, the quantities ^'^^T^^ which appear on the left side 
of Eq. (7), arc the Christoffcl symbols associated with the induced metric 
in the leaves of foliation defined above. 



3.1 Case 1: p^l 

For the class of warped geometries given by Eq. (4) and (5), we can easily 

see that for p = 1 we have ^'^''r^^ = and (5)r^4 = y''" (%^)- 

In this case, the 4-dimensional part of the geodesic equations reduces to 
the form 

(Px^^,{4)rtid^d^_nf>dx'^dy 
dX^ + -0 dX d\ ~ ^ dX dx- ^ 
where a prime denotes differentiation with respect to y. Similarly, the geodesic 
equation for the fifth coordinate y in this warped product space becomes 

d'^y w 9f , dx" dx'^ „ 

^+,/'e^^/..,^^=0. (11) 

The above two geodesic equations (10) and (11) are identical to those 
in Ref. [17]. This happens due to the special nature of the warping function 
chosen by us. We like to point out that the type of warping function in our 
case is different from that used in [17], as they considered a warping function 
which depends on the extra coordinates only. So, we can conclude that 
the geodesic equations do not depend on the explicit form of the warping 
function for this particular choice when p = 1. The result can be stated in 
the form of the following proposition: 

Proposition! : The geodesic equations for the 4-dimensional spacetime 
and for the fifth dimension remains the same as in the case with a warp 

factor which depends only on the extra coordinate, even after the inclusion 
of the additional linear time-dependence of the warping function, provided 
9yy — l- 



A study of the 5-dimensional motion of particles with non-zero rest mass 
near the hypersurface can be done by considering the 5-dimensional timelike 
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geodesies {9ab^§^^jx' ^ which the above equation can be easily 

decoupled from the 4-dimensional spacetime coordinates as in [17], yielding 



0. (12) 



The motion of photons in 5-dimension can similarly be analyzed by 
considering the null geodesies, for which ((?ab %^ ~ ^) ^^'^ ^® from 
Eq. (11) 

'''+/'f'lV = 0. (13) 



dA2 ■' VrfA, 

Equations (12) and (13) are second order ordinary differential equations, 
which can be solved if we know the nature of the warping function /. 



3.1.1 Motion of particles in four- dimensional spacetime To determine the 
nature of motion of particles in the four-dimensional world, we need to find 
the constraints on the warping function that will determine whether the 
geodesies on the bulk manifold may coincide with those on the hypersur- 
faee or not. Considering the expression for the extrinsic curvature of the 
hypersurface in this case, we find that the result is identical to those in Ref. 
[18]. 



3.1.2 Motion of particles in the fifth dimension A qualitative analysis of 
the motion in fifth dimension can be done without actually solving Eqs. (12) 
and (13) by defining q = ^ and investigating the autonomous dynamical 
systems [37] 



^^Pii^y) (15) 

with P{q, y) = — ./'(e + 9^), where e = 1 for timelike geodesies and e = for 
null geodesies. The equilibrium points of the system of equations (14) and 
(15) are given by ^ = and = 0. Knowledge of these points along with 
their stability properties can provide a lot of information about the behavior 
allowed by this dynamical system. For the warping function considered by 
us, the analysis is very similar to the one done by Daliia et. al. [17]. 

Therefore, we find that the results obtained in [17] and [18] are true even 
if we include a linear time-dependence of the warping function, provided 
9yy — l- 
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3.2 Case 2: p=p(t,y) 

The geodesic motion of the test particle in the bulk spacetime is still de- 
scribed by the equations (6). Let us consider a bulk metric of the form: 

dS^ = e^f^^'y^ [dt^ - btdr^ - btr^de^ - btr^sin{efd(P'^) - p{t, y)dy'^ (16) 

where b is another constant and / is given by Eqn. (5) as f{t,y) = at + 
lln{cosh{cy)). On substituting the explicit expression for the bulk metric, 
the geodesic equations reduce to the following five second order difi^erential 
equations: 



dtY^{2at + l)b 



d'^r {2at + 1) dt dr 
d>? 2t dXdX'"^ 





sm 



+ sm - 



\dXj 



dX 



0, 



= 0, 
(17) 

(18) 



rf26» d9 r {2at + 1) dt , ldr\ 
d)^^dX\ 2f 



= 0, 



(19) 



d^ d4 
dJ^^dX 



{2at + 1) dt 1 dr cos 9 d9 
2t dX r dX sinO dX 



= 0, 



(20) 



rfA2 2p\ dt dX dX dy\dx' ^ 



(21) 



Let us now analyse the equations of 5-dimensional motion of both pho- 
tons as well as the particles of non-zero rest mass, near the hypersurfacc. For 
that we consider the equations of the 5-dimensional geodesies to be given 

by 



dzA dz^ 



where, e = 1 for timelike geodesies and e = for null geodesies. Using the 
expression for the bulk metric, we get from (22), 



dr 
dX 



Substituting from (23) in (17), we get 



p2/ 



(23) 
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+ 



{2at + 1) 
2t 




0, (24) 



Once again, a qualitative analysis of the motion of particles can be done 
without actually solving the Eqs. (17) to (21) by defining a dynamical sys- 
tem in the following way: 



and 



U = 



V 



dt 
dX' 



dy 



which, when substituted in (24) and (21) yields the equations 

1 \ 1 



dU 



U'[2a 



2t 



o2f 



0, 



and 



dX 2p\ dt 



0. 



(25) 



(26) 



(27) 



(28) 



The set of equations (25) to (28) are used to analyse the evolution of the 
above dynamical system. The trajectories of this dynamical system can be 
described in a 2-dimensional phase space in terms of the variables U and V. 
The system of equations (25) to (28) represent a real non-linear dynamical 
system [38], [39] of the type 



dt 
dX 
dy 
dX 
dU 

Ix 

dV 

Ix 



u, 

V, 

P{U,V,t,y), 

QiU,V,t,y), 



(29) 



where 



P{U,V,t,y) = -U^ (2a + 



QiU,V,t,y) 



1 



1 



1 f dp , ^^29p 



+ ^{e+pV^} 



dy 



with e = 1, for timelike and null geodesies respectively. Both P and Q have 
continuous first partial derivatives for all (U,V). The evolution of the system 
in the (U,V) phase plane is given by the solution of the set (29). It is evident 
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that the solution of the system for a given value of the parameter A is no 
longer unique, thereby, giving rise to very complicated trajectories. In such 
cases, the critical point of the system must be a zero of the equations for 
all times for which the system is defined. Since no linear terms are present 
in the set (29), there will generally be several critical points of the system. 
We know that the zeros of the system correspond to the fixed points of the 
phase trajectories. For a nonautonomous system, generally each point of 
the phase space is intersected by many distinct trajectories. Therefore the 
collection of phase trajectories as well as the analysis for the identification 
of the fixed points becomes extremely complicated [40] . The dynamics of the 
phase trajectories is then strongly controlled by the underlying topology of 
the manifold under consideration. Such an analysis is beyond the scope of 
the present paper. One simple method is to make suitable transformations to 
reduce the non-autonomous systems to autonomous ones by appending t to 
the depending variables. However, such a transformation has the drawback 
that the non-autonomous problem frequently looses its original algebraic 
structure. Hence it is not considered here. 



4 Phase trajectories and Critical Points of the system for the 
case p=p(t,y) 



In spite of the fact that we are dealing with a nonautonomous system, we 
can still interpret the nature of the trajectories with the help of a simple 
analysis. The equation 

^ - ^ (30) 
dU~ 

at some given value of y specifies the phase trajectory of the system in the 
(U,V) phase plane, provided P 7^ at this point. Generally P ^ except 
at the point U = V = ioi null geodesies (which, as we shall find later, 
represents the critical point of the system in such a case). Substituting the 
explicit forms of P and Q in (30) and integrating, we arrive at the result 



-pIt-^''+2t\ 



where K is the arbitrary constant of integration. Setting this constant to 
zero and simplifying, we can recast eqn. (31) into the form 

AU'^ + BUV + CV'^ + D = 0, (32) 

where A, B, C, and D are functions of coordinates, specifically, t and y. The 
above equation is useful to understand the phase trajectories even when no 
fixed points are found. 

For null geodesies, D = and (32) reduces to 

AU^ + BUV + CV^ = 0. (33) 
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For a given value of A, B and C, the phase trajectories are represented by 
a pair of straight Hnes intersecting at (0,0). The slope of the lines change 
for different choices of A, B and C. It is evident that for null geodesies, 
P = = Q at U = V = Q, which therefore represents the fixed point of the 
trajectories in this case. 

However, for timelike geodesies, P ^ even when U = V = 0. In that 
case, eqn. (32) represents the phase trajectory of the system even if no fixed 
point is found. In fact, we shall find from our subsequent analysis, that we 
cannot locate any fixed point for timelike geodesies, except for a special 
case. In the following figures, we have illustrated a few possible trajectories 
for some arbitrarily chosen values of the functions A, B, C, and D in case 
of timelike geodesies. The Fig. 1 represents the plot for the sample equation 
W^-GUV+SV^ = 2, Fig. 2 represents the equation -U"^ ~6UV +8V'^ = -2, 
Fig. 3 represents the equation U"^ — 6UV — SV"^ = —2.5 and Fig. 4 represents 
the equation 4U^ — 6UV + 8V^ = 100. We can interpret the different figures 
to correspond to different values of y and hence they lie on different (U,V) 
phase planes. A particular trajectory lies on a given (U,V) phase plane with 
y = constant. 



Figl . eps 



Fig. 1 Plot for the sample equation U"^ - QUV + SV'^ 
The trajectory appears to avoid the point (0,0). 



2 for timelike geodesies. 



Initial examination of the above figures, indicates a few possibilities. 
Namely, the elliptic trajectory represents a closed orbit with U = y = as 
the center and thus an example of neutral aquilihrium,. Such an orbit repre- 
sents a system whose total energy is conserved. In general, the trajectories 
appear to avoid the point U = V = 0. 

Thus when the specific warp factor is coupled with a lapse function, 
which depends on both time and extra coordinates, the physical system is 
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Fig2 . eps 



Fig. 2 Plot for the sample equation —U^—6UV+8V^ = —2 for timelike geodesies. 
This trajectory also appears to move away from (0,0). 



Fig3.eps 



Fig. 3 Plot for the sample equation U'^-QUV-W'^ 
This trajectory also behaves like the earlier ones. 



-2.5 for timelike geodesies. 



modeled by a nonautonomous dynamical system. This is different from the 
results obtained for p = 1, as well as from those obtained by Dahia ([17] and 
[18]), who dealt with an autonomous system. We can therefore say that, the 
physical system in this case is a more general one. 



We now proceed to examine the critical points of the system analytically. 
Let {Uo,Vo) represent the critical point. 
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Fig4 . eps 



Fig. 4 Plot for the sample equation 4U^—6UV+8V^ = 100 for timelike geodesies. 
This is a special case indicating a closed orbit, that signifies neutral equilibrium. 



Null geodesies: In this case, the two conditions P = and Q = lead 
us to 



2a 



1 \ , PVS 



and 



From (34), we have 



2t e2/ 



Vq = Uq 

Substituting from (36) in (35), we get 









' 2t) 







2W^ + W'^P 



dt 



dy. 



= 0, 



(34) 



(35) 



(36) 



(37) 



where, W = and is not zero in general. Thus, in general, 

(37) will be valid if we have Uq = Q and hence, Vq = 0. Thus (0,0) will be 
the critical point for null geodesies. 



Timelike geodesies: In this case, the condition P = = Q lead us to the 
result 
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1 (opY 
P \9y J 



(38) 



w^^(l)'-4(i?) 



2 ' 



and 





(39) 



P 



Owing to our assumption regarding the nature of p, neither of Uq or Vq 
can be zero in general, nor can be imaginary unless 



Thus, the critical point cannot be determined unless the exact form of 
the function p{t, y) is known. In fact, there will not be any fixed point unless 
the condition P = = Q is maintained for all times for which the dynamical 
system is defined. 

5 Cosmological Interpretations 

The cosmology of the 4-dimensional hypersurfaces is determined by com- 
puting the Einstein tensor directly from the metric corresponding to the 
5-diniensional spacetime. We have used the GRTensor package [41] for our 
calculations. For the type of line element given by (4) and (5), we choose 
the particular case. 



The 5-dimensional metric in general represents an observable universe 



with isotropic pressure given by Pr = pe = P4, = j [ e(2,f+2*tf(coffe(c«)))t2 



Both ordinary matter in the bulk as well as the induced matter on the 
hypersurface exhibit isotropic pressure given by the above expression. The 
geometric properties of the bulk spacetime remains unchanged for all values 
of c including zero. 

Proposition2: It is possible to have matter residing on the correspond- 
ing 4-dimensionaI hypersurface in the form of pure radiation. The geometric 
properties of the bulk spacetime remains unchanged for all values of c in- 
cluding zero. 




(40) 




There are two important points to note. First, for a = 0, the warping 
function reduces to the type given by Gremm [42]. In that case we have 
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Gab = Rab, and the Ricciscalar is identically zero. The 5-dimcnsional Ricci 
tensor is then identical to the Ricci tensor of the corresponding 4-metric 
and the matter on the hypersurface is in the form of pure radiation with 
G( = and = Gg = G'^ = Secondly, the Einstein tensor for the 

5-metric given by (16) is identical to the 5-metric obtained by the trans- 
formation ln{cosh{cy)) 1 and hence the kind of 4-dimensional geometry 
observed here can also be embedded in other 5-metrics obtained by such 
transformations. 

It is interesting to note that the observable universe still represents mat- 
ter with isotropic pressure even in the case where the lapse function is a 
function of both time and extra coordinates. 

6 Conclusions 

The consideration of the confinement of the geodesies of particles and the 
cosmology associated with the corresponding 4-dimensional hypersurfaces 
is an important program in the braneworld scenario. In this paper we 
have done it for RS-type braneworlds with the bulk in the form of a five- 
dimensional warped product space-time. The warping function is an expo- 
nential function of both time and extra coordinates and the lapse function 
is either a constant or a function of both time and extra coordinates. In 
such a case, we have obtained a mechanism for the confinement of geodesies 
on co-dimension one hypersurfaces. 

Considering that confinement is purely due to classical gravitational ef- 
fects and geometric in nature, we have obtained a description of the geodesic 
motions for these warped product spaces by using the natural splitting that 
occurs between the motion in the extra dimension and the motion in the 
four- dimensional hypersurfaces. Further, this splitting helps us to use the 
phase space analysis to determine the nature of the geodesic motions in the 
neighbourhood of hypersurfaces. It is observed that the geodesic equations 
for tlie 4-dimensional spacetime and for the fifth dimension do not depend 
on the time dependence of the warping function provided the lapse hmc- 
tion is a constant. Further, in the general case when the lapse function is 
a function of both time and extra coordinates, using the dynamical sys- 
tem analysis we have determined the nature of the trajectories and have 
examined the critical points on a given hypersurface. 

Prom the point of view of cosmology, the observable imiverse in general 
represents matter with isotropic pressure and which, in a particular case 
will be in the form of pure radiation. 

For future work it is interesting to study whether the time dependence 
of the warp factor is related to the time evolution of the dark energy. 
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